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Abstract. We define an integrable hamiltonian system of Toda type associ- 
ated with the real Lie algebra so(p,q). As usual there exists a periodic and a 
non-periodic version. We construct, using the root space, two Lax pair repre- 
sentations and the associated Poisson tensors. We prove Liouville integrability 
and examine the multi-hamiltonian structure. The system is a projection of a 
canonical A n type Toda lattice via a Flaschka type transformation. It is also 
obtained via a complex change of variables from the classical Toda lattice. 



1. Introduction 

The Toda lattice, is a Hamiltonian system with an exponential potential: 



N N-l 



N , N—i 

H(q 1 ,...,q N , p 1 ,...,p N ) 

i=l ~ i=l 

It describes a system of N particles on a line, connected by exponential springs. 
The Toda lattice is a well-known integrable system with soliton solutions. Under 
Flaschka's change of variable [5] 



(1) 



a = i e 4(«-«+0 
2 



b, = -~Pi , 



the equations become 



(2) 



<k {bi+i - h) 
Hal -all) 



These equations can be written in Lax pair form L = [B,L], where L is the 
Jacobi matrix 
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and B is the skew-symmetric part of L in the Lie algebra decomposition lower 
triangular plus skew-symmetric. Due to the Lax pair, it follows that the functions 
Hi = j tr L 1 are constants of motion. Moreover, they are in involution with respect 
to a Poisson structure, associated to this Lie algebra decomposition. 
We note that 



N 1 
F!=5> = --(pi 



-P2 



Pn) 
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corresponds to the total momentum and 



N JV-1 



H 2 = H(q u ...,q N ,pi,... ,p N ) = \^2,b 2 + ^2 a i 

»=l »=i 

is the Hamiltonian. 

Consider R 2N with coordinates (g±, . . . , Qn,Pi, ■ ■ ■ ,Pn), the standard symplectic 
bracket { , }o and the Flaschka' transformation F : R 2N — > R 2W_1 defined by 

F : (q 1 ,...,q N ,p 1 ,...,p N ) -t (01, . . . , a N -i, h, . . . , b N ) . 

There exists a bracket on R 27V_1 which satisfies 

{f,9}°F = {foF,goF} . 

It is a bracket which (up to a constant multiple) is given by 

{a l: b l+1 } = ^ ; 

all other brackets are zero. Hi = b\ + 62 + • ■ • + is the only Casimir. The 
Hamiltonian in this bracket is Hi = h tr L 2 . We also have involution of invariants, 
{Hi, Hj} — 0. We denote this bracket by ttl. 

The quadratic Toda bracket, Ii 2 , due to Adler [I], can be used to define a bi- 
hamiltonian formulation of the system. It is a Poisson bracket in which the Hamil- 
tonian vector field generated by Hi is the same as the Hamiltonian vector field 
generated by Hi with respect to the tti bracket. The defining relations are 



{ai, at+i} 
{ai,bi} 



1 



^ {ai,b. l+ i} = aib i+ i 



{b h b i+ i} = 2 a?; 

all other brackets are zero. This bracket has det L as Casimir and Hi = tr L is the 
Hamiltonian. Furthermore, H2 is compatible with 7Ti. We also have the Lenard 
relations 

(6) n 2 Vif; - mVH l+ i . 

The classical Toda lattice was generalized to produce a Toda type system for 
each simple Lie algebra. The finite, non-periodic Toda lattice corresponds to a 
root system of type Ag. This generalization is due to Bogoyavlensky [2]. These 
systems were studied extensively in |10) where the solution of the system was con- 
nected intimately with the representation theory of simple Lie groups. We call 
these systems the Bogoyavlensky- Toda lattices. They can be described as follows. 

Let q be any simple Lie algebra over the complex numbers. One chooses a Cartan 
subalgebra t) of g, and a basis n of simple roots for the root system A of f) in g. 
The corresponding set of positive roots is denoted by A + . To each positive root 
a one can associate a triple (X a , X- a , H a ) of vectors in g which generate a Lie 
subalgebra isomorphic to s^C). The set (x a , X-a) Q <=A+ U (^a)aen is basis of g, 
called a root basis. To these data one associates the Lax equation L = [B, L], where 
L and B are defined as follows: 



L = hh ai + a,i(x ai + x-gj, 

i=l i=l 

i 

B = ^ai(x at - x_q,J. 
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The affine space of all elements L of q of the above form is the phase space of 
the Bogoyavlensky-Toda lattice, associated to q. The functions which yield the 
integrability of the system are the yld-invariant functions on g, restricted to M. 

In this paper we imitate this procedure for the case of a real Lie algebra of type 
so(p,q). The resulting system is Liouville integrable as well. It turns out that this 
system can be realized as a projection of a classical type Toda system where some 
of the terms in the potential are switched from positive to negative. We illustrate 
with a small example: The potential of the Toda lattice is exponential 

jV-l 

V(q u ...,q N ) = Y / z q *~ q * +1 ■ 
i=l 

For example, if N = 5 

We may modify the potential by changing some of the plus sings to minuses the 
system remains integrable. There are 16 total possibilities for the possible ± signs 
in front of the exponentials. For example one may consider the potential 

V(qi,q 2 ,q 3 ,qi,q 5 ) = e« + - e «»-* +e «"* . 

By projecting this system using Flaschka's transformation we obtain a new sys- 
tem in (a, b) variables and we will give a Lax pair for this system. It turns out 
that the Lax pair is associated with the real Lie algebra so(6,5). More generally 
we define an integrable hamiltonian system associated with the real Lie algebra 
so (2m, 2n + 1). As usual there exists a periodic and a non-periodic version. 

We construct, using the root space, a Lax pair representation and the asso- 
ciated Poisson tensors. We prove Liouville integrability and examine the multi- 
hamiltonian structure. The system is bi-hamiltonian with a linear and a quadratic 
Poisson bracket as in the case of the classical Toda lattice. The quadratic Poisson 
bracket is similar to the Adler Toda bracket with some differences on the signs. In 
fact one can obtain this new bracket via a simple complex change of variables. The 
system is a projection of a canonical A n type Toda lattice (with some sign changes 
on the potential) via a Flaschka type transformation. It turns out that this system 
is equivalent via a complex change of variables with the classical non-periodic Toda 
lattice. As we mentioned, there are a total of 2 N such systems but in this paper 
we just consider one such potential for simplicity. Note that the number of systems 
is related to the construction of Tomei [15] of the Toda manifold, i.e. the space of 
real, symmetric, tridiagonal nxn matrices with fixed eigenvalues. 

We begin by giving a description of the basis of so(2m, 2n + 1), its Cartan 
subalgebra, its roots and its root space. We end-up with a new set of polynomial 
equations in the variables (a, b). One can write the equations in Lax pair form 
(L(t), B(t)), which can be described in terms of the root system. We also compute 
the Casimirs associated with both Poisson brackets and prove the involution of the 
invariants. Finally we show that this system is equivalent via a simple complex 
change of variables to the classical Toda system. Using this approach we give a 
second Lax pair involving complex coefficients and a different set of integrals in 
involution. 

2. Bl-HAMILTONIAN SYSTEMS AND MASTER SYMMETRIES 

A Poisson bracket on the space C°° (M) of smooth functions on a smooth man- 
ifold M is a skew-symmetric, bilinear map, 

{•, •} : C°°(M) x C°°(M) -> C°°(M), 
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that verifies the Jacobi identity and is a biderivation. Thus, (C°°(M), {•, •}) has 
the structure of a Lie algebra. A Poisson structure on a smooth manifold M is a 
Lie algebra structure on C°°(M) whose the bracket {•, •} : C°°(M) x C°°(M) -> 
C°°{M) verifies the Leibniz's rule: 

{/, 9h} = {/, g}h + g{f, h}, V /, g, h e C°° (M). 

The Poisson bracket {•, •} gives rise to a contravariant antisymmetric tensor field 7r 
of order 2 such that ir(df,dg) = {/, g}, for f,g£ C°°(M). The bivector 7r is called 
a Poisson tensor and the manifold (M, tt) a Poisson manifold. 

A Hamiltonian system is specified by a Poisson bracket { , } together with a 
C°° function H called the Hamiltonian. The vector field has the form 

X = ndH . 

A bi-Hamiltonian system is defined by specifying two Hamiltonian functions Hi, 
H2 satisfying 

(7) X = 7Ti dH 2 = tt 2 dff a . 

If 7Ti is symplectic, i.e., the structure matrix tt± is invertible, we call the Poisson 
pair non-degenerate. If we assume a non-degenerate pair we make the following 
definition: The recursion operator associated with a non-degenerate pair is the 
(1, 1) -tensor 1Z defined by 

(8) n = 7T 2 TT- 1 . 

We have the following result due to Magri |12) : 

Theorem 1. Suppose that we have a non- degenerate bi-Hamiltonian system on 
a manifold M , whose first cohomology group is trivial. Then, there exists a hier- 
archy of mutually commuting functions Hi,H 2 , ■ ■ ■ , all in involution with respect 
to both brackets. If we denote by \% the Hamiltonian vector field generated by Hi 
with respect to the initial bracket %i then the Xi generate mutually commuting bi- 
Hamiltonian flows, satisfying the Lenard recursion relations 

(9) Xi+j = n V-ffj , 

where iti — TV^ni are the higher order Poisson tensors. 

The definition and basic properties of master symmetries can be found in Fuchssteiner 
[S]. Consider a differential equation on a manifold M defined by a vector field X- 
A vector field Z is a symmetry of the equation if 

(10) [z, x ]=o. 

If Z is time dependent, then a more general condition is 

37 

(11) — + [Z,x} = 0. 
A vector field Z is called a master symmetry if 

(12) [[Z,x},x}=0, 
but 

(13) [Z,X]*0- 

Master symmetries were first introduced by Fokas and Fuchssteiner in [7] in 
connection with the Benjamin-Ono Equation. 

Suppose that we have a bi-Hamiltonian system defined by the Poisson tensors 
7Ti, 7T2 and the Hamiltonians Hi, H 2 . Assume that w± is symplectic. We define the 
recursion operator 1Z — 7T27rf , the higher flows 

(14) Xi=K i ~ 1 xi, 



SO(P,Q) TODA SYSTEMS 



and the higher order Poisson tensors 

For a non-degenerate bi-Hamiltonian system, master symmetries can be gener- 
ated using a method due to W. Oevel [14] . 

Theorem 2 (Oevel). Suppose that Xq is a conformal symmetry for both ~k\, 1T2 
and Hi, i.e., for some scalars X, fi, and v we have 

Then the vector fields 

Xi = TZ 1 Xq 
are master symmetries and we have 

(a) C Xt H 3 =(y+(j-l + t)( p - \))H i+j , 

(b) C Xi TTj = (p + (j-i- 2)(/i - X))n+j , 

(c) [X i ,X j } = (n-\)(j-i)X i+j . 

3. The real Lie algebra so(p, q) 

In this section we construct a Toda type system associated with the real alge- 
bra so(p,q). The construction is similar in spirit to the one of Bogoyavlensky for 
complex simple Lie algebras. 

We start with the real Lie algebra so(p, q) where p — 2m and q = 2n + 1 for 
some positive integers m,n. 

To simplify the notation we let N = m + n. Recall that by definition (see [IT]) 

S0(p, q) = { X G Ql(p + q, M) | X l I v . q + I p>q X = } 
where X 1 is the transpose of X and 



I P 
-I q 



We represent the elements of so(p, q) by real (p + q)-by-(p + q) matrices of the 
form with a and d skew symmetric, of size p-by-p and q-by-q respectively. 

A basis for so(p, q) consists of the matrices e%$ for l^i<j^p + q where: 

(a) For 1 ^ i < j ^ p and p+l^i<j^p + q 

( 1 if (k,l) = (i,j) 

(15) the fc/-th entry of = < — 1 if (k, I) = 

{ if (k,l)^(i,j) 

(b) For 1 ^ i ^ p and p+l^j^p + q 

( 1 if (k,l) = (i,j) 

(16) the fcZ-th entry of = < 1 if (k,l) = (j,i) 

{ if (k,l)^(i,j) 

A basis of the Cartan subalgebra f) consists of N block-diagonal matrices of size 

/ 1 N 

2N + 1, whose first N blocks are of size 2-by-2, with i-th block J = I ^ _ 

and whose last block is of size 1-by-l. Note that we can express the elements of 
f) in Kronecker product form if we drop their last row and column which consist 
entirely of O's. Hence, let di^i = diag(0, . . . , 1, . . . , 0) be of size N-by-N with 1 in 
the i-th position and let ® denote the Kronecker product of matrices. Then the 
i-th basis element of F) is given by di,i ® J and is equal to en-\,u with the last row 
and column removed. 
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We denote the basis elements of f) by h ak where the ak, are the simple roots 
defined as follows. For 1 ^ k ^ N — 1: 

!(— l) fe i j = k, 
(~l) fe i j = k+ 1, 
otherwise. 

Alternatively the simple roots correspond to vectors of size N as follows: 

ai = (— i, — i, 0, . . . ,0), 
a 2 = (0,i,i,0,...,0), 
(18) "3 = (0,0, — i, — i, 0, ...,0), 

a N -i = (0,0,...,(-l) N -H,(-l) N -H). 

The root vectors corresponding to the roots ak, —ak, for 1 ^ k ^ N — 1, are given 
by: 

(IQ) Xa k = 1/2 (e2fe-i,2fe+i — e2fe,2fc+2 + (-l) fc (e2fc-i,2fc+2 + e2fe,2fc+i)i) 

x- ak = 1/2 (e 2 fe-i, 2 fe+i - e 2 fe,2fc+2 + (-l) fe+1 (e2fc-i,2fc+2 + e2fe,2fc+i)i) • 
We can express the root space in Kronecker product form as follows. Firstly, let 

(20) J k = [ { _l )k . M f ^ for 1 < fc < N - 1. 

Secondly, for 1 < fc < N — 1, let -Xj, be an N-by-N matrix defined as follows: 
For k ^ m: 

( 1/2 if (i,j) = (k,k + l) 

(21) the ij-th entry of X k = < -1/2 if = {k + l, k) 

[ otherwise. 

For k = m: 

1/2 if = (m,m+ 1) 

(22) the ij-th entry of X m = { 1/2 if (i, j) = (m+l, m) 



otherwise. 



Finally, for 1 ^ k ^ N — 1, we have: 



One can define the system in Lax pair form (L{t),B{t)) in 5, using the root 
system as follows: 

JV JV-l 

L(t) = £(-l)'6j (t)h aj + (-iya 3 (t)(x aj + x- aj ) , 
(23) ^V-x 

B W = -iEH)^WK-t- a ,)' 

j'=i 

As usual /i aj is an element of a fixed Cartan subalgebra and x a . is a root vector 
corresponding to the simple root aj. It is straightforward to verify that the Lax 
equation L(t) = \B{t),L{t)] is consistent and gives the following equations: 

= l,2,...,iV-l 



(24) 



hi 


= ^ (b l+1 - b{), 


k 


= 2(o?-a?_i), 


h 


= 2a? 


b m 


= -2 ^ + a 2 m _ 


b m +i 


= 2(a? n+1 +a^) 


b N 


= -24-1 ■ 
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Before giving an explicit matrix form of the Lax pair we remark that the last 
row and column of the matrices L and B consist entirely of zeros. For the sake of 
simplicity we will omit the last row and column of L and B since we still have a 
consistent Lax pair. 

The Lax pair (|23p . can be written explicitly (after omitting the last row and 
column) using the Kronecker product as follows: 



(25) 



where 



(26) 



L = Li® 



Li 







-1 

1 





-b 2 



1 

-1 



\ 



(-l) N+1 b N J 



(27) 



/ -a x 
a\ a 2 
-o 2 '•• 



Ln = 





; o (-i^-Wi 

V ••• •••0 (-1)^0^-! o 

Note that the upper diagonal of L 2 is 

(-01, 2 , ■ • ■ , (-l)*Oi, . . . , (-l^V-l), 

whereas the lower diagonal is 

(oi, -02, • ■ • , (-l) m o m _i, (-l) m a m , (-l) m+2 a m+1 , (-l) N a N ^). 
That is 



(28) the ij-th entry of L 2 



(-l)X ifj = i + l, 

(-iy +1 aj i£j = i-l,j^m, 

(-l) j aj iij = i-l,j = m, 

otherwise. 



The matrix B is defined as follows: 



(29) 



B = Bi 



1 

1 
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where 



( ai 
— ai a 2 
-a 2 



(30) 



B x = 



\ 



a m •• 

: •• •• •• ajv_i 

\ ••• ••• -a N -i / 

Note that the upper diagonal of B\ is (ai, a 2 , . . . , a/v-i), whereas the lower diagonal 
is 

( — 11, • • • j — Ctm-l; a m, -dm+l) • ■ • j — OjV-l)- 

It turns out that the trace of odd powers of L arc all zero. 
Lemma 1. 

H 2t +i =trL 2l+1 =0 Vi . 
Proof. Recall that by definition 

ao(p, g) = { X g fll(p + <7,R) | + 7 Pi<z X = } 

where X* is the transpose of X and 







It follows that X 
i.e., 



-I p ^X t I p ^ q . Let /„(x) be the characteristic polynomial of X, 
/ n (A)=det (A7-X) . 



Then 



/„(-A) = det(-AJ - X) = (-l)"det(A7 + X) 

= (-l)"det(A/ - I-^XH^ = (-l)"det(7- g 1 7 P;(! 
= (-l)"det(J-^(AJ - = (-l) n dct(A7 - X 1 ) 



-ly 



Since the Lax matrix is of odd dimension, n is odd. Dividing by A (which cor- 
responds to eigenvalue zero) we obtain an even polynomial. Therefore, for the 
Lax matrix L if A is an eigenvalue of L so is —A. This implies that characteristic 
polynomial of L is an even polynomial. It also shows that tr L k = if k is odd. □ 

As we will see shortly the function Hi = b\ + b 2 + • • • + frjv-i (which is not the 
trace of L) will turn out to be a Casimir. Therefore we need N — 1 independent 
functions in involution is order to establish integrability. 

The following functions are constants of motion: 

(-iy 



H. 



2i 



Ai 



txL z \ for i = 1,...N- 1. 



4. SYMPLECTIC REALIZATION 

In R 2N with coordinates (q, p), q = (q\, . . . , gjv), and p = {pi,... ,Pn) we choose 
the following Hamiltonian function: 



(31) 



N 1 

%-p) = Eo? 




aii-gi+i 



-9m+l 
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which differs from the classical, non-periodic Toda lattice by a minus sign at position 

TO. 

Using Flaschka's transformation F : (q, p) — > (a, b) we obtain equations (|24|) . 
Recall that the (a, 6) variables are defined by: 

( o 2] ^ = ^(ft-am), i = l,2,...,N-l, 

1 1 h = -\pi, i = l,2,...,N. 

Let the mapping F : R 2Ar — > R 2JV_1 , where 

F : (qi,...,qN,Pi,---,PN) -» (ai, . . . , ajv-i, &i, • • • , 6jv) , 

denote the Flaschka transformation given by equation ([1]). Recall that the stan- 
dard symplectic bracket on R 2JV is mapped onto the linear bracket tti in (a, b) 
coordinates. The linear bracket has only one Casimir Hi = bi + b 2 + ■ • • + b^- 
We can take as Hamiltonian the function H 2 = — \ tr 1? , where L is given by 

We define a quadratic bracket which we call tt 2 with defining relations: 
{ai,bi} = -dibi , 



1,.. 


..,N- 


1, 


1,.. 


..,N- 


1, 


1,.. 


..,N- 


1, 


1,.. 


. . , TO — 


l,m + l,...,N - 1, 



(33) {a i7 6 i+ i} = a^j+i , 
{6i, 6*+i} = 2 a 2 , 

{&m>&m+l} = — 2a m ; 

all other brackets are zero. Note that this bracket is different from the Adler bracket 
defined in ([3]). We remark that the mapping of R 2Ar_1 — > R 27V_1 given by 

(a 1 ) -> (a 1 ,...,a m _ 1 ,m m ,a m+ i,...,ajv_i,&i,6 2 ,...,6„) 

sends the Adler bracket II2 onto the bracket 7T2 . Since n 2 is known to be Poisson this 
implies that tt 2 is also Poisson. We simply think of the two brackets as complex 
brackets connected by a complex Poisson map and therefore the Jacobi identity 
holds. Since the coefficients of the defining functions of 7r 2 are real it is also Poisson 
as a real Poisson manifold. 

This bracket also has one Casimir det L and H± = b\ + . . . bjy is the Hamiltonian. 

Furthermore, tt 2 is compatible with tti. We also have 

(34) TT^H-2 = lT 2 VHl . 

Similarly we define another Poisson bracket, tt^, which is a modification of the 
cubic Toda Poisson bracket. The defining relations for W3 are 



(35) 



{at, a i+1 } 


— aiai + ibi+i , 


i = l,.. 


.,N- 


2, 








{a,i,bi} 


= - ai b1 - of , 


i=l,.. 


. , TO — 


1, to 4 


-1,... 


,N- 


1, 


{a,i,bi} 


= -a t b 2 + of , 


i = to, 












{ai,b i+ i} 


= ai b 2 l+1 + a? , 


1 = 1,.. 


. , TO — 


1, m -\ 


-1,... 


,N- 


1, 


{di,b i+ i} 


= aib 2 l+1 - af , 


i = to, 












{a>i, h+2} 


2 

— a i a i+l 7 


i = l,.. 


. , TO — 


2, TO, . 


..,JV 


-2, 




{ai, b i+2 } 


2 

— ~ a i a i+l ! 


i = to — 


1, 










{a l+ i,bi} 


= -afa l+1 , 


i = l,.. 


. , TO — 


1, m -) 


-1,... 


,iV- 


2, 


{a l+ i,bi} 


2 


i = to, 












{bi,b i+ i} 


= 2a 2 (b % + b z+1 ) , 
= ~2a 2 {h + b l+l ) , 


i = l,.. 


. , TO — 


l,mH 


-1,... 


,N- 


1, 


{bi,b i+1 } 


i = to; 













all other brackets are zero. 

The bracket 1T3 is compatible with both tti and 7T 2 and satisfy Lenard type 
relations like 

TTldH^ = TT^dHl . 
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We note that 

N 



(36) H 1 =^rb i = --{p 1 +p 2 + ---+p N ) 

i=l 

corresponds to the total momentum (Hi ^ tr L) and 



N ( N-l 

\^+[ E 

i=l \i=l,i^ 



N-l 

2 I „2 



(37) H 2 =H(qi,...,q N ,p!,...,p N ) = ^ E &i + E a * I ~~ a 
is the Hamiltonian. 

We now define Toda hierarchies for the Toda lattice in (q,p) variables. We follow 
reference [I]. 

Let Ji be the standard symplectic bracket with Poisson matrix 

J> 
-I 0, 



J i = 



where I is the N x N identity matrix. We use J\ — 4 j\. With this convention the 
bracket J\ is mapped precisely onto the bracket ir\ under the Flaschka transforma- 
tion (JTJ) . We define J2 to be the tensor 

A 

-B C , 



■h 



where A is the skew-symmetric matrix defined by a%j = 1 = — aji for i < j, B 
is the diagonal matrix (— p±, —P2, • • ■ , ~Pn) arid C is the skew-symmetric matrix 
whose non-zero terms are c^+i = — Cj+i^ = e qi ~ qi+1 for i — 1, 2, . . . , m — 1, m + 
l,...,iV- 1 and c m , m+ i = -c m+ i jm = - e 9"»-?m+i_ We define J 2 = 2J 2 . With 
this convention the bracket J2 is mapped precisely onto the bracket 712 under the 
Flaschka transformation. 

It is easy to see that we have a bi-Hamiltonian pair. We define 

h\ = -2(pi +p 2 H \-pn) , 

and I12 to be the Hamiltonian: 



N ( N-l 

E^+ E 

i= 1 yz__l,'i^:r. 



Under Flaschka's transformation (JXJ) , /ix is mapped onto 4(&i + 62 + • ■ ■ + &w) = 
47?i and /12 is mapped onto — trL 2 = AH2- Using the relationship (|34p . we obtain, 
after multiplication by 4, the following pair: 

J 1 Vh 2 = J 2 V/ii . 

We define the recursion operator as follows: 

n = J2J1 1 • 

The matrix form of 1Z is quite simple: 

< 38 > K = Kc ~B 

In (q,p) coordinates, the symbol Xi i s a shorthand for x^. It is generated as 
usual by 

Xi = n i - 1 xi . 

In a similar fashion we obtain the higher order Poisson tensors 
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We finally define the conformal symmetry 

N r) N 1 

Zo = g (A ,_ 2j + 1) _ + g B _. 

It is straightforward to verify that 

£z Ji = — Ji j 
Cz J2 = . 

In addition, 

Z (hi) = hi 
Z (h 2 ) = 2h 2 . 

Consequently, Zq is a conformal symmetry for Ji, J 2 and h\. The constants 
appearing in Oevel's Theorem are A = — 1, /j, = and v = 1. Therefore, we end-up 
with the following deformation relations: 

[Zi,hj] = (i + j)h l+j 
CziJj = (j — i — 2)Ji+j 
[Zi,Zj] = (j - i)Z i+j . 

We note that Magri's theorem implies the Lenard type relations 

(39) Jidhj — Jjdhi . 

It is straightforward to verify that the symplectic bracket J\ is mapped via 
Flaschka's transformation F onto the linear Poisson bracket n\ and J 2 is mapped 
onto the quadratic bracket it 2 . Similarly hi, h 2 correspond via F to H\ and H 2 
respectively. The master symmetries Zi are mapped via F to produce a sequence 
of master symmetries Xi in the (a, b) phase space. 



5. Integrability 

In this section we prove the integrability of the so(p, q) Toda system. We note 
that the tensors J\ and J 2 project under Flaschka's transformation to the linear 
and quadratic brackets ~k\ and ~k 2 respectively. As in the case of the classical Toda 
we have the bi-hamiltonian formulation. 

(40) 7Tl Vff 2 = 7T 2 Viil . 

Of course this bi-hamiltonian formulation is quite different than the Lenard relations 
(|39p . The Poisson tensors tt\ and tt 2 are no longer symplectic and as we know there 
exists no recursion operator in the (a, b) space. It is well-known (see [5], [4], [13] ) 
that the recursion operator 7Z cannot be reduced. It is also well-known that the 
vector fields Zi are projectable and give rise to a sequence of vector fields X\ , X 2l . . . 
in the (a, b) space [5J. To prove the involution of the invariants we will only use the 
master symmetry X 2 and the reduced Hamiltonians H 2 i which are the constants 
of motion for our system. We will make use of the Lenard relations 

(41) TridH 2l = -n 3 dH 2l - 2 Vi . 

Theorem 3. The functions H 2 i are in involution with respect to the tt\ Poisson 
bracket. 

Proof. Using the Lenard relations (|4T|) it is straightforward to prove the involution 
of the functions H 2i . The Poisson bracket we use {, } is the linear Poisson bracket 

7Tl. 
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Since H 2 is the Hamiltonian we have {H 2 ,H 2 j} = Vj since H 2 j is clearly a 
constant of motion. Now we calculate 



{H 2i ,H 2 j} — 



< dH 2l , -K\dH.<xi > 

- < dH 2j , TTidH 2l > 

- < dH 2j , n 3 dH 2i - 2 > 

< dH 2 i-2, K3dH 2 j > 

< dH 2 i-2, nidH 2j+2 > 
{H21-2, H 2 j +2 } 



{H 2 , H 2i+2 j^ 2 } — 



□ 



It is not difficult to see that the functions H 2 i, i = 1,2, . . . , N — 1 are functionally 
independent on an open dense set. Using the original approach of Henon [S] it is 
enough to assume a\ = a 2 = ■ ■ ■ = a^-i = 0. Then the traces of L 21 are symmetric 
functions of the hi. In fact taking into account the form of the Lax matrix (1251 ) it 
is easy to see that 



H 2i = 



+ b 21 



N I 



at a = 0. Therefore the Jacobian matrix has the form 



( 5 1 


b 2 


bl 


bl 


b\ 


b\ 



\b\ N ~ z b™-y 

Selecting the top left corner N — 1 x TV — 1 minor we obtain a Vandermonde 
matrix whose determinant is different than zero on the set bi ^ bj, for i =/= j. 

We can prove integrability in a different way, by finding a different Lax pair and 
a new set of independent functions in involution. Dchnc the following Lax pair 
(A, M) where 



A ? 

b% 



,2N-3 














\ 






ai 


b 2 


a 2 




(42) 


M = 





a-2 


b 3 














io m : 












ia m ■■ aAr_i 






U 









i.e. we replace a m with ia m where i 2 = —1. We take A to be the skew-symmetric 
part of M. One verifies easily that the Lax pair is equivalent to the equations ([24]) . 
Therefore the functions 

Ii = T trM\ i = l,2,...,iV 
i 

are first integrals of the system (j2~4"|) . As we have seen I\ = R\ = 61 + b 2 + ■ ■ ■ + ojv 
is the Casimir of the Poisson bracket 7Ti . Therefore we may use the set of integrals 
{I 2 , 13, . . . In} as a different set of functions to be used to prove integrability. Note 
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that Hn is equal to In ■ To prove that this set is in involution one can use the 
methods of [5] . Recall that the first master symmetry of the classical Toda lattice 



is 






(43) 


^ a»n — 

n—l 


N a 

n—l 


where 






(44) 


OL n = -na n b n - 


f (n + 2)o„6„+i 


(45) 


(3 n = (2n + 3)a 2 + 


(1 - 2n)<£_i + bl 



Replacing a m by ia m in the first equation does not alter the equation. On the 
other hand f5 m now changes to 

p m = -(2m + 3)a 2 „ + (2m - l)^ + b 2 m . 

With this change we obtain master symmetry X\ which has the following prop- 
erties. 

X 1 (H j ) = (j + l)H j+1 

C Xl TT 3 = (j - 3)7T J + l j = 1, 2, 3 . 

One now establishes the integrability of system (l24l) following [3l |4] . The proof 
of involution of invariants is exactly the same as in [U Proposition 3] . 

6. Example so(6, 5) 

We illustrate in detail the results with a specific example. 

Using the root structure of so (6, 5) we define the following Lax pair. 

The matrix L is given by 



L = 



( 


-bi 


-ai 




















o \ 


bi 








ai 




















ai 








b 2 


a 2 




















-ai 


-b 2 




























-a 2 








-b 3 


-a 3 




















2 


b 3 








a 3 




















-a 3 








bi 


a 4 




















f'3 


-bi 








— 04 




















— 0,4 








-h 


V o 




















a 4 


h 


o ) 



or equivalently in Kronecker product form 



L 



(0 


-ai 










ai 





a 2 











-a 2 





-a 3 











-a 3 


a 4 










—0,4 


°/ 


(bl 








0\ 







-b 2 
















63 













- 


-64 












6 5 y 







-f 



-f 





The matrix i3 is the following: 
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B = 



I 











ai 




















\ 










ai 






























-ai 











a 2 


















-a\ 











a 2 






























-ai 











«3 


















-ai 











«3 






























a 3 











a 4 


















a 3 











CL 4 






























— <24 










\ 




















— Cl 4 











/ 



or equivalently 





( ° 


ai 








°\ 




-CLl 













B = 





-a 2 





03 


















a 4 




V 








— (24 


0/ 



The Lax pair is equivalent to the following equations of motion: 



hit) 
hit) 
hit) 
hit) 

The constants of motion are Hn — 4i 
to establish intcgrability since Hi 
For example: 



2 fll 2 

2a 2 2 - 2ai 2 
-2a 3 2 - 2a 2 2 
2a 4 2 + 2a 3 2 
-2a 4 2 



01 (t) 

02 (*) 

03 (*) 
a' 4 (t) 



(62 
ih 
(64 
(65 



61) ai 

62) 02 

6 3 ) a 3 

64) 04 



(-1)' 



tr L 2j i = 1,2,3,4 which are sufficient 



Sj=i i s a Casimir for the 7Ti bracket. 



H 4 



= h+b 2 + b 3 + b 4 + 65 



1 



= ai 2 + a 2 2 - a 3 2 + a 4 2 + -(6 2 + 6 2 2 + h 2 + h 2 + b 2 ) 

= ai 2 (fei 2 + &1&2 + &2 2 ) + a2 2 (&2 2 + Ms + h 2 ) - a 3 2 (h 2 + hh + h 2 ) 
+a 4 2 ib 4 2 + b 4 b 5 + b 5 2 ) + a 2 a 2 2 - a 2 2 a 3 2 - a 3 2 a 4 2 
+l/2( fll 4 + a 2 4 + a 3 4 + a 4 4 ) + l/4(6i 4 + b 2 4 + b 3 4 + b 4 4 + b 5 4 ) 

Note that H 2 , H 4 , H§,H% form a set of independent integrals of motion in invo- 
lution demonstrating the integrability of the system. 

We obtain another set of integrals in involution in the following way: Define the 
Lax pair (M, A) where 



(h 

a 1 



M = 



ai 

h 





a 2 

a 2 h 
ia 3 







ia 3 
h 



0\ 





0,4 



a 4 b 5 J 



and 
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( 


ai 








o\ 


-ai 





a-2 











-ai 





ia 3 











-ia 3 







V o 








— (24 





Define: Ji = Hi, - (/,•_! ) 
We have = 2J7 2 . Define 

^2 = --X"i(.ffi) = • 

Then define 

J 3 = ^i(ff 2 ) , 

and 

J 4 = ^i(J 3 )=fl4 • 

Finally let 

J 8 = ^i(4) • 
5 

Then the set {/ 2 — Hi, h,l4 = H4,I 5 } is another set of independent integrals 
in involution. 
Of course 

L = - tr M j . 

j 

For example 

1 5 

In addition we have {Ij,Hk} = for j, k = 1, 2, 3, 4. Therefore the functions i?6 
and i7 8 must be functions of the 7, . For example 

77s — H^-\-H-\Ik-\-HoHa 77?77H — U?Iv — 77i 77 2 7 3 HfHo — 77?H — 7*?H — H 

720 1 2 1 6 1 24 1 6 2 2 3 4 

The equations of motion can be also obtained using the Lie-Poisson bracket 7Ti 



7Tj = 



/ 














-on 


ai 











\ 





















ai 






























-03 


03 






























— a a 


a 4 






a-i 






























-ai 


a-2 






























-a 2 


03 






























-03 


04 



















V 











—04 

















/ 



and the Hamiltonian i7 2 - 

There is one Casimir for m: 

Hi = h+b 2 + --- + b 5 . 
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The quadratic Poisson bracket is defined by: 



7T 2 



u 


a±a,2 
2 


u 


u 


— a\Oi 


ai0 2 


u 


u 


n 
u 


2 





2 








-a 2 6 2 


a2^3 











2 





a3a4 
2 








-a 3 6 3 


a 3 b 4 











a 3 a 4 
2 














— a 4 b 4 


a4& 5 


aih 














2 fll 2 











~aib 2 


a 2 b 2 








-2 ai 2 





2a 2 2 











-a2&3 


a3^3 








-2a 2 2 





-2a 3 2 











-a 3 6 4 


0464 








2a 3 2 





2a 4 2 











-a 4 6 5 











-2a 4 2 






V 

There is one Casimir for 7r 2 : 

V det L = —bib 2 b 3 b 4 b5 + bib 2 b 3 a 2 — b\b 2 b^a 3 + 616465(12 — b\a\a\ 
+b 3 b 4 b^a\ — 6 3 a 2 a 2 + b^a\a\ = — detM 

The formulas for the master symmetries X\ and X 2 are the following: 

/ §01(561-62) \ 
\a 2 (3 6 2 + 6 3 ) 
\ a 3 (3 6 4 + 6 3 ) 
- \ a 4 (64 — 5 65) 
Xi = -2ai 2 + 6i 2 

4 ai 2 +6 2 2 
2a 2 2 -2a 3 2 + 6 3 2 
4a 4 2 + 6 4 2 
V -2a 4 2 + 6 5 2 J 

\ a\ (2 6i6 2 + 6461 + 6561 - 6 3 6 2 — 6 4 6 2 — 6562 + 6 3 6i + 5 6i 2 — 6 2 2 + 2 ai 2 ) 
\ a 2 (3 6 2 2 - 6i6 2 + 2 6 3 6 2 + 6 4 6 2 + 6 5 6 2 + 4 ai 2 + 2 a 2 2 - 2 a 3 2 + 6 3 2 + 6 3 6i - 6 4 6 3 - 6 5 6 3 ) 
6361 - 6 3 6 2 + 2 6 4 6 3 + 6 5 6 3 + 4 a 4 2 + 3 6 4 2 + 6 4 6i + 6 4 6 2 - b 5 b 4 + 2 a 2 2 - 2 a 3 2 + 6 3 2 ) 
\ a 4 (6 5 6 3 + 65^2 + 6561 - 6 4 6 3 - 6 4 6 2 - 6 4 6i + 2 6564 + 2 a 4 2 - 6 4 2 + 5 6 5 2 ) 
-2 6 2 ai 2 - 6iai 2 - ai 2 6 3 - ai 2 6 4 - ai 2 6 5 + 6i 3 
4 6iai 2 + 5 6 2 ai 2 + «i 2 6 3 + a 1 2 b 4 + ai 2 b 5 + a 2 2 bi + b 2 a 2 2 - a 2 2 b 4 - a 2 2 b 5 + b 2 3 
2 b 2 a 2 2 - a 2 2 h + 3 b 3 a 2 2 + a 2 2 b 4 + a 2 2 b 5 - 2 b 4 a 3 2 - a 3 2 b x - a 3 2 b 2 - 3 b 3 a 3 2 + a 3 2 b 5 + b 3 3 
a 3 2 bi + a 3 2 b 2 - b 4 a 3 2 - a 3 2 b 5 + 4 b 5 a 4 2 + a 2 b x + a 2 b 2 + a 2 b 3 + 5 b 4 a 2 + b 4 3 
\ -2 b 4 a 2 - a 2 b x - a 2 b 2 - a 2 b 3 - b 5 a 4 2 + b 5 3 

7. Generalizations 

We conclude with two possible generalizations of the systems considered. 

There is also a periodic version of the system which is also integrable. To obtain 
it we make some simple modifications in the various definitions. In this section 
we indicate briefly the modifications needed to obtain the periodic version of the 
system. 



X 2 



2 

2-«3 
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The last root, a^, is given by: 
(46) 



i j = 1, 



a N (h aj )={ {-l) N i j = N, 

otherwise. 

Alternatively the simple roots correspond to vectors of size N as follows: 

o-l = (-L -i,0, ...,0), 



(47) 



Q' 2 = (0, i, i, 0, . . . ,0), 

a 3 = (0,0,-i,-i,0,...,0), 



a N = (i,0,...,0,(-l) jv i). 

Again we obtain a Lax equation L(t) — [B(t), L(t)] which is consistent and gives 
the following equations for 1 ^ i ^ N: 

hi = ai(b i+ i ~bi), for i ^ N, 
a N = a N (h - b N ), 



(48) 



bi 
h 

K, 



2(a?-a?_ 1 ), for i <£{l,m,m+ 1,N}, 
2 (a? + a%) 
-2(<& + a£,_i) 



— 2 (a m _)_i 



= -2(a 2 N + a 2 N _ x ) . 
In the periodic case in (q, p) coordinates we use the Hamiltonian 



(49) 



N I N-l 

tf(q,p)=£-P? + E e9 ^ 9,+1 



-<? m +i 



which is slightly different from the classical, periodic Toda lattice. (The minus sign 
at position m and the last term.) In order to obtain equations (14"51) . in the variables 
(a, b) , we extend the Flaschka type transformation given in Q32p by adding a new 
variable a at. 

(50) a N = i e ^ 9N ~ 9l) . 

The Lax pair L = [B, L] can be written using the Kronecker product as follows: 

(51) L = Li®fi ~ { f)+ L 2®(^ _°i) +i 3®^ (_!)JV+i 
where L\ and L 2 are as in (|2"B")) and (|2"7)l respectively and 

/ (-l) w a^\ 

: '•• 



(52) 



La = 







(53) 



\(~l) N a N 
The matrix B is defined as follows: 

1 

1 



/ 



B = B 1 



Bo 



1 

{-l) N 
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where B\ is as in (|3T)1) and 

/ 



(54) 



B 2 = 



a N \ 




\(-l) N a N 







0/ 



There is also a different Lax pair (M, A) 







/6i 


a i 





iaw \ 






ai 


b 2 


2 




(55) 


M = 





a, 2 
















ia m : 












ia m '■ ajv-i 






\iOAT 






otv-i &jv / 



i.e. we replace a m with ia m . We define 





(56) 



.4 



/ ai 

— Ol 02 

-a 2 



\iajv 



-iajv\ 



ia m : 

OjV-1 

-Oat-i / 



One verifies easily that the Lax pair is equivalent to the equations (148)) . 
Note that the Hamiltonian has one more term in (a, b) coordinates i.e. 

^ n I N-l 

(57) H 2 =H(qi,...,q N ,pi,...,p N ) = -^2b?+ \ < 



' N 



Of 

(58) 



The Poisson bracket wi for the periodic case is defined as in (fj| with the addition 

{ojv, h} = a N . 



There are two Casimirs Hi = b± + 62 + • ■ • + fe/v and 0102 ■ ■ • a^r. 
There is a quadratic Toda bracket tt2 with defining relations as in (|33p with the 
addition of 



(59) 



{ojv, oi} = — \a,Nai , 

{a*N,bN} = —a-NbN , 

{a N ,bi} = a N bi , 

{b N M} = 



-2a% ; 



This bracket also has two Casimirs: det L and 0102 ■ ■ ■ aj\f. 
Hi = b\ + . . . 6jy is the Hamiltonian. 
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Furthermore, -k 2 is compatible with -k\. We also have the Lenard relation 

(60) 7TiViJ 2 = 7r 2 Viii . 

We can generalize these systems in a different way by choosing different signs in 
the potential terms. Tomei in [T5] studies the topology of the set of real, symmetric, 
tridiagonal n x n matrices with fixed eigenvalues 

Ai > A2 > • • • > Ajv ■ 

This set consists of 2 Ar ~ 1 components depending on the signs (±) of the variables 
ax, 0,2, ■ ■ • , We may construct a Toda system corresponding to each component 
and any two of them are isomorphic via a complex transformation. In this paper we 
have studied in detail one such example but in general the number of such systems 
is 2 Ar ~ 1 . It is also possible to produce two Lax pairs for each system. The first one 
is obtained from the Flaschka Lax pair by replacing some of the aj by iaj. The 
second Lax pair which is obtained in so(p, q) is also easily obtained by modifying 
some of the signs of (|23|) . We illustrate with an example for N = 3. 

Example 1. In this example we display all possible Toda systems in the case of 
N — 3. We give a Lax pair and the equations of motion for each system. For the 
last system we indicate also the second Lax pair which is obtained by an obvious 
change of signs in i25\) . 



Example 2 



Li = 


bi ax 
ax b 2 a 2 
\0 a 2 63 


of motion: 




b\(t) 


= 2 (a?) 


b 2 (t) 


= 2(ai-c 


h(t) 


= 2(-a|) 


2. 




L 2 = ( 


fbx ax 
ax b 2 ia 2 
V ia 2 63 


of motion: 




b\{t) 


= 2 (a?) 


Wt) 


= 2(-a|- 


b 3 (t) 


= 2 (a|) 





f 


a 1 







-ax 





«2 




\ 


-02 






dx(t) 
d 2 {t) 



(b 2 - b x )ax 
(63 - b 2 )a 2 





( 


ax 





B 2 = | 


-a x 





i«2 




\ 


— ia 2 






a\(t) 

02 (t) 



(6 2 - bx)ax 
(63 - b 2 )a 2 



Example 3. 



La 




Bo 






iax 










«2 





-02 






Equations of motion: 
b\(t) 
b 2 (t) 

h(t) 



2 (-a 

2(al 

2 (-a 



01 (i) = 

02(*) = 



(62 - b x )ax 
(63 - 6 2 )a 2 
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Example 4. 











iai 


•°1 











iai 


■°» 
















iai 




ici2 




-iai 





lfl2 
















,o 


ia 2 


63/ 











-ia 2 / 










Equat 


ions of motion. 
































b\(t) 




2 (-a 


1) 




01 (t) 






(62- 


- 6i)ai 


















2 (-a 


i + a?) 




02(t) 






(63- 


- ^2)02 


















2(a|) 
























Alternative Lax pair for the last example: 




















( ° 


-h 


-Oi 








\ 






/ 








ai 








\ 




h 








Oi 

















01 













L> 4 : 


-ai 








&2 


a 2 





#4 









ai 








a-2 







ai 


-62 





- 


-0,2 






01 








a-2 




















- 


-b 3 














02 












{ ° 








-a 2 


h 


° J 






V 





a 2 











/ 
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